Journal of Nonlinear Mathematical Physics 1998, V.5, N 3, [28J4293 



Article 



> 



X 



A Solvable Many-Body Problem in the Plane 

F. CALOGERO 

* Dipartimento di Fisica, Universita di Roma "La Sapienza" , 00153 Roma, Ital^ 
Istituto Nazionale di Fisica Nucleare, Sezione di Roma, Raly 



OO 

Q\ • Received May 1, 1998 

ON ' 



Abstract 



A solvable many-body problem in the plane is exhibited. It is characterized by 
rotation-invariant Newtonian ("acceleration equal force") equations of motion, fea- 
turing one-body ("external") and pair ( "interparticle" ) forces. The former depend 
quadratically on the velocity, and nonlinearly on the coordinate, of the moving parti- 
cle. The latter depend linearly on the coordinate of the moving particle, and linearly 
respectively nonlinearly on the velocity respectively the coordinate of the other par- 
' ticle. The model contains 2n 2 arbitrary coupling constants, n being the number of 

particles. The behaviour of the solutions is outlined; special cases in which the motion 
C ■ ' is confined (multiply periodic), or even completely periodic, are identified. 

O ' 
OO 

1 Introduction 

Recently several solvable and/or integrable many-body problems in the plane have been 
^ introduced [l]-[4]. (We use here the heuristic - imprecise but useful - distinction among 

solvable models, for which the solutions can be obtained in relatively explicit form, and in- 
tegrable models, which possess a sufficiently large number of independent globally defined 
constants of the motion). In this paper we report one more such model, which seems 
worthy of special notice because of its simplicity and generality (2n 2 arbitrary coupling 
constants, see below). Indeed its equations of motion read as follows: 



2fj{fj -Tj) -fjjfj -fj) 
' 3 



(1-1) 

mST(R m„ ^ h ^k{r k ■ rj) + rj{r k ■ r k ) - r k {r k ■ rj) 



r 2 

k=i k 



Notation. Superimposed arrows identify 2-vectors, for which we use the notation 

f=(x,y), zAf=(—y,x), r 2 = x 2 +y 2 . (1-2) 
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(The second formula corresponds to the 3-dimensional notation f= (x,y,0), z = (1,0,0), 
with the symbol A denoting the standard 3-dimensional vector product). The index j 
always takes integer values from 1 to n, the total number of interacting particles. The 2n 2 
"coupling constants" (5jk and 7^ are arbitrary. Dots denote of course differentiation with 
respect to the time t. 

The equations of motion (1.1) are of Newtonian type ("acceleration equal force"), with 
one- and two-body forces. The one-body ("external") force is quadratic in the particle 
velocity fj, and depends nonlinear ly on the particle coordinate fj. The two-body ("in- 
terparticle" ) force depends linearly (only) on the coordinate of the moving particle, and 
linearly respectively nonlinearly on the velocity respectively the coordinate of the other 
particle. The model is rotation-invariant, but not translation-invariant (see however be- 
low for a translation-invariant version). It is invariant under the scale transformation 
fj — ► cfj, with c an arbitrary constant. 

This model is solvable, as explained in the following Section. In the subsequent Section 3 
a brief discussion is given of the particle motions it entails, and the following generalized 
version of (1.1) is also introduced: 

Z /, ~ x 4 2fj(fj ■ fj) — fjifj ■ fj) 



r) 

(1.3a) 



k=i Tk 
Pjk(t) = [Pjk cos(ut) - ~f jk sin(wt)] exp(At), 

(1.3b) 

ljk{t) = [7jfcCos(wt) + (3 jk $m(ut)\ exp(At). 

For A = ijj = 0, this generalized model reduces to the previous one; for A = 0, w / 
(and of course real), all its solutions are completely periodic, with period T = 2tt/uj (see 
Section 3). In Section 3 we also indicate how to manufacture analogous solvable models, 
which are also translation-invariant. 

Let us finally emphasize that many-body models in the plane generally exhibit a much 
more interesting behavior than models on the line; this is confirmed by the results reported 
below, as well as by those displayed in previous papers [l]-[4]. 

2 The solution 

In this Section we indicated how to solve the system (1.1). 

Let the n complex quantities fj(t) evolve in time according to the linear firstorder 
evolution equations 

n 

fj = ^2 a jkfk, (2.1) 
k=l 

whose general solution reads 

n 

/i(t) = X>fexp(a fc t), (2.2) 
k 
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where the n quantities a k are the n eigenvalues of the (n x n)-matrix A with elements ctj k , 
and the n n- vectors <3?W, with components <p( k \ are the corresponding eigenvectors: 



E(m) (m) / Q \ 

a jkfk = a m¥) (2.3) 

k=l 



Here we are implicitly assuming, for simplicity, that the matrix A is diagonalizable, and 
that all its eigenvalues a k are distinct; the diligent reader will have no difficulty in detailing 
the modifications in the treatment which are required if these assumptions do not hold. 
Introduce now the n quantities zj(t) via the position 

f j = z j /z j , (2.4) 

which clearly entails 

f j = z j /z j -z]/z], (2.5) 
hence, via (2.1) and (2.4), 

n 

Zj = Zj/zj + ^2a jk ZjZ k /z k . (2.6) 
k=l 



On the other hand, from (2.4) and (2.2), one gets 



Zj(t) = Zj(0)exp 



X^ fc) [ ex P(°fci) - 1] /a k 
.k=i 



(2.7) 



while the normalization of the n n- vectors <3?W, with components (pf \ is fixed by the 
conditions 



f>f =i,(0)/^(0), (2.8) 
k=i 

also implied by (2.2) and (2.4). 

Hence (2.6) is explicitly solvable, for any choice of the n 2 quantities aj k , and 2n initial 
data Zj(0) and ij(0). And this is true as well if all these quantities are complex, so that 
we can put 

z j = x j +iy j , (2.9) 

ajk = Pjk + iljk, (2-10) 

of course now with Xj, yj, (3j k and jj k all real. 

But it is now a trivial exercise to verify that, via these positions, the complex evolu- 
tion equations (2.6) coincide with the 2- dimensional rotation-invariant real equations of 
motion (1.1), whose solvability is thereby demonstrated. 
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3 Discussion of the motions, and other models 

In this Section we briefly discuss the motions yielded by the many-body problem (1.1), 
using the findings of the preceding Section(see in particular (2.7)); and we then indicate 
how to obtain other solvable models in the plane, including the many-body problem (1.3) 
and translation-invariant generalizations of the many-body problems mentioned above. 

Clearly the most important element that characterizes the behavior of the system (1.1) 
are the n eigenvalues of the matrix A. If all these eigenvalues have negative real parts, 
from any initial conditions the system will tend to a standstill configuration; all particle 
velocities vanish exponentially as t — > oo. Such an outcome may also occur, but only for 
special initial conditions, if only some of these eigenvalues have negative real parts. A 
necessary and sufficient condition for the system to possess periodic trajectories is that 
at least one of these n eigenvalues be imaginary. If all these eigenvalues are imaginary, 
the motion remains confined, and is in fact multiply periodic, for any initial conditions. 
If moreover all these eigenvalues are rational multiples of the same imaginary quantity, 
then for any initial conditions the system behaves completely periodically. If one of these 
eigenvalues vanishes, the system (2.6) possesses the similarity solutions 

Zj(t) = Zj(0)exp(rit), (3.1a) 

with rj = A + iuj an arbitrary (complex) constant. The corresponding similarity solution 
of (1.1) reads 

fj(t) = exp(At) [cos(wt) + sin(otf)zA] r}(0). (3.1b) 

If some of the n eigenvalues a k have positive real parts, the generic solution may feature 
some (from to n) particles which escape, doubly exponentially fast, to infinity, while the 
others converge, doubly exponentially fast, to the origin. 

Let us now derive the system (1.3a) with (1.3b). To this end let us start from the 
system (2.6), but with the quantities Zj(t) formally replaced by, say, Cj( T ) ( an d accordingly 
with dots replaced by, say, primes, the latter indicating derivatives with respect to r). We 
then introduce new quantities Zj(t) via the position Zj(t) = so that (2.6) gets 

replaced by 

n 

Zj = Zj (f/f) + zj/zj + f a jk z jZk /z k . (3.2) 

k=l 

Now we make again the transition from the complex dependent variables Zj(t) to the real 
2-vectors fj(t) (as above, see (2.9) and (1.2)). It is then easily seen that the new particle 
coordinates fj(t) satisfy precisely (1.3a) with (1.3b), provided 

r = {exp{(\ + iuj)t] - l}/(A + iw). (3.3) 

The statement made after (1.3b) is thereby proven. 

Finally, let us manufacture a generalized version of (1.1), which is translation-invariant. 
To this end we introduce a model with 2n particles of two types, whose coordinates we 
denote respectively by rj + \t) and rj \t), and we set 

r 3 (t) = fj + \t) - r}~\t), Rjit) = f} + \t) + r}~\t), (3.4) 
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with fj(t) satisfying the (solvable) equations of motion (1.1), and Rj(t) satisfying the 
(clearly, also solvable; see below) equations of motion 

Rj = (Aj + n j zA)R j . (3.5) 

Here Xj and Clj are In arbitrary (of course real) coupling constants. 

We thereby see that the new "particle coordinates" r^\t) evolve according to the 
translation-invariant equations of motion 



^■ ±) = 2{( A i + n ^ A )4-± 



r) 



(3.6) 



fc=i rfc 



of course with (3.4). 

As for (3.5), a convenient way to exhibit its solution is via the formula 

Zj{t) = Zj(0) + Zj(0) {exp [(Aj + iilj)t]} /(A, + ify), (3.7) 

which clearly provides the general solution of (3.5) via the identification of the complex 
dependent variables Zj = Xj + iYj, with the real 2- vectors Rj = (Xj, Yj), since this 
identification entails that, to (3.5), there corresponds 

Zj = (Aj + iQj)Zj. (3.8) 

Clearly the same trick may be directly applied to (1.3a) with (1.3b), yielding a tran- 
slation-invariant generalization of these equations of motion. An alternative possibility, 
generally yielding a different model, is to apply the change of independent variable (3.3) 
to the complex scalar version of (3.6), and subsequently to change the dependent variables 
(from scalar complex numbers to real 2-vectors). We do not exhibit the corresponding 
results, which the interested reader may immediately write down. Nor do we report an 
analysis of the actual behavior of the solutions of (3.6) with (3.4), or of the new models 
just mentioned, since the above results are so transparent not to require elaborations. But 
we end by reemphasizing the reachness of the motions entailed by these formulas, which 
is now enhanced by the additional freedom connected with the choice of the 2n constants 
Aj and fl,. 

References 

[1] Calogero F., A Solvable n-Body Problem in the Plane. I, J. Math. Phys., 1996, V.37, 1735-1759. 

[2] Calogero F., Three Solvable Many-Body Problems in the Plane, Acta Applicandae Mathematicae, 
1998, V.51, 93-111. 

[3] Calogero F., Tricks of the Trade: Relating and Deriving Solvable and Integrable Dynamical Systems, 
in: Proceedings of the Workshop on Calogero-Moser-Sutherland Models, Montreal, 10-15 March, 
1997 (in press). 



[4] Calogero F., Integrable and Solvable Many-Body Problems in the Plane via Complexification, J. 
Math. Phys. (in press). 



